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The onset of intruder ground states in Na isotopes is investigated by comparing experimental 
data and shell-model calculations. This onset is one of the consequences of the disappearance of the 
N=20 magic structure, and the Na isotopes are shown to play a special role in clarifying the change 
of this magic structure. Both the electromagnetic moments and the energy levels clearly indicate 
an onset of ground state intruder configurations at neutron number A*' = 19 already, which arises 
only with a narrow = 20 shell gap in Na isotopes resulting from the spin-isospin dependence of 
the nucleon-nucleon interaction (as compared to a wider gap in stable nuclei like ''"Ca). It is shown 
why the previous report based on the mass led to a wrong conclusion. 
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I. INTRODUCTION 

Amongst the most intriguing and unique features in ex- 
otic nuclei are rather significant changes from the conven- 
tional magic structure. As a result of them, the ground 
state of a nucleus with N or Z close to a conventional 
magic number is not necessarily spherical, and can be 
strongly deformed. Its fingerprint was first identified 
from extra binding energies of "^^■'^^Na 0], whose origin 
was regarded, consistently with a Hartree-Fock calcula- 
tion y , as the dominance of strongly deformed intruder 
components in the ground state over the normal compo- 
nents. Here, normal (intruder) states imply the states 
comprised of shell-model configurations without (with) 
Iplh, 2p2h or higher excited configurations across the 
iV = 20 shell gap. Later, more direct experimental evi- 
dence of the strong deformation was found for ^^Mg from 
the low excitation energy of the 2f state 0, Q and the 
large B{E2; 0+ -> 2+) Q value. Thus, the disappear- 
ance of the A^=20 magic structure has been established 
in some = 20 isotones including the recent case for 
•^"^Ne 6]. It still remains, however, an open question as 
to where the ground state changes from a normal- to an 
intruder-dominant configuration in the chain of isotopes, 
and the question as to what mechanism drives this disap- 
pearance remains. The present paper aims at presenting 
the resolution of these questions, as exemplified in the 
structure of Na isotopes. 

For Na isotopes, one may expect that the onset of 
the intruder-dominance of the ground state lies right at 
A = 20, from the comparison of the experimental mass 
to a shell-model result within the sd shell with the USD 



interaction 0. The USD interaction has been the most 
frequently used interaction in the sd shell, and we shall 
refer to shell model calculations with this interaction in 
the sd shell as USD model or calculation, hereafter. A 
similar picture about the onset is assumed in the so-called 
"island of inversion" model |^ |^ , where the lowest nor- 
mal and the lowest intruder states are confronted with- 
out mixing between them. Although the mass (or the 
separation energy) often provides us with helpful infor- 
mation on shell structure, studies from different angles 
are needed before one draws definite conclusions, as we 
shall demonstrate. The first part of the present paper is 
focused upon re-examination on the dominant configura- 
tion of the ground state of Na isotopes. We perform a 
large-scale shell-model calculation using the Monte Carlo 
shell model (MCSM) 10], which is bricfiy described in 
Sec. m In Sec. IIIII and IIVI respectively, the electro- 
magnetic moments and the energy levels are discussed, 
and from such discussions the transition point from the 
normal- to intruder-dominant ground state is identified 
in the chain of Na isotopes. In Sec. El the second part 
of the present paper, we discuss the mechanism of the 
disappearance of the magic structure, focusing upon the 
(effective) A^=20 gap between the sd and pf shells and 
emphasizing the special importance of the nucleus '^'^Na 
on this issue from a somewhat general viewpoint. We 
finally summarize the present study in Sec. IVII 
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II. OUTLINE OF THE SHELL MODEL 
CALCULATION 

The model space and the effective interaction used in 
the present study are the same as those of our previous 
studies [ill ll^ : the valence shell consists of the full sd- 
shell orbits and two lower p/-shell orbits. The effective 
interaction is called hereafter SDPF-M for the sake of 
clarification from other interactions. The SDPF-M inter- 
action was introduced in Ref. TT] in 1999, by combining 
the USD interaction (7J for the sd shell, the Kuo-Brown 
interaction for the p.f shell, and a modified Millener- 
Kurath interaction for the cross shell. On top of 
this, a small but important modification was made for its 
monopole part |^l| as we shall add some remarks later. 
A unique feature of the SDPF-M interaction is that the 
neutron shell structure, defined by the effective single- 
particle energy (ESPE), chang function of the pro- 
ton number, more significantly than in previous models, 
for instance, the "island of inversion" Here, the 
ESPE includes mean effects from other valence nucleons 
on top of the usual single-particle energies with respect 
to the given inert core (i.e., closed shell). Therefore, the 
ESPE depends on shell-model interactions between va- 
lence nucleons. The present varying shell structure can 
be explained by the shell evolution mechanism of Ref. 
[TH in terms of the spin-isospin property of the effective 
nucleon-nucleon (NN) interaction. The strong T = 
monopole attraction between the 0^3/2 and the Od^/2 en- 
larges the = 20 gap, as protons occupy 0^5/2. In- 
versely, this effect diminishes towards Z — 8, ending up 
with a rather narrow A = 20 gap and a wider A = 16 
gap. This shell evolution leads us to the oxygen drip line 
at A = 16 [ll|p,[ll[ll as a resuh of emerging A = 16 
magic number l2(J| . The monopole part of the SDPF-M 
interaction was modified from that of the USD interac- 
tion so as to reproduce the oxygen drip line pl| , while the 
resultant monopole part is closer to the G-matrix result 
as emphasized in jl5j . 

In the A = 20 region, as we shall illustrate, the Na 
isotopes give indispensable information on this shell evo- 
lution: (i) specific Na isotopes provide us with clues 
of a narrow A = 20 shell gap, (ii) with odd Z, their 
ground-state properties can be directly examined by non- 
vanishing electromagnetic moments, and (iii) many ex- 
perimental data have been recently accumulated about 
the mass moment ,2^221: 7-ray spectrum and tran- 
sition by the Coulomb excitation |24L |25| . etc. Thus, we 
carry out shell- model studies on Na isotopes from A=16 
to 20. 

Since the dimension of the Hamiltonian matrix be- 
comes prohibitively large with the present problems, we 
perform a shell-model calculation by the MCSM based 
on the quantum Monte Carlo diagonalization (QMCD) 
method whose development has been described in |2a . l27l 
[28, 29] . In the present MCSM calculation, we adopt the 
so-called J -compressed bases p^ . i.e., bases generated 
and adopted by monitoring the energy with the full angu- 



lar momentum projection. The feasibility of the MCSM 
calculation for odd-A nuclei has been demonstrated in 
[l^. This method works very well for odd-odd nuclei as 
well. 

In the present calculation, the E2 matrix elements 
are calculated with the effective charges (ep,e„) = 
(1.3e,0.5e) which are the same as those used in the 
USD model 0- It has been confirmed that the MCSM 
with these effective charges excellently reproduce the 
B{E2;0t 2+) values of even-even nuclei from stable 
to unstable nuclei for the effective Ml oper- 

ator. Brown and Wildenthal took an empirically opti- 
mum one within the USD model [s^]- They found that 
the free-nucleon g factors give no obviously deviating 
magnetic moments but more quantitative agreement can 
be attained with the empirically optimum operator: for 
A = 28, gs is quenched by a factor 0.85, and gf = 1.127, 
3;" = -0.089, gP = 0.041 and g^ = -0.35 are used where 
I, s, and p are the orbital angular momentum, the in- 
trinsic spin, and V87r[F(^^(r) (g) s]'-'^' operators, respec- 
tively. Recently, Honma et al. have presented in [sj 
that the spin g factor does not have to be much quenched 
in the full p/-shell model space using their newly devel- 
oped interaction |^ [s^ • Based on these extensive shell- 
model studies, the g factors are adopted, in the present 
MCSM calculations, so as to be rather close to the above- 
mentioned ones. Namely the spin part is quenched by a 
factor 0.9, and the other g factors are shifted from the 
free-nucleon values by SgiilY) = 0.15 and SgpilV) = 0.5 
where g(IV) denotes the isovector g factor defined by 
5(IV) = (gP-5")/2. 



III. ELECTROMAGNETIC MOMENTS 

The shell-model calculation described in the above sec- 
tion is carried out for Na isotopes from A = 16 to 20. We 
first compare, in Fig.^ the electric quadrupole moments 
and the magnetic dipole moments between the MCSM 
with the SDPF-M interaction and experimental data. As 
a reference, results from the USD model are presented, 
also. For the A = 16 and 17 isotopes, the experimental 
moments [23, 113, H^] are well reproduced by both the 
shell-model calculations, reflecting the dominance of the 
sc?-shell configurations in their ground states (see Fig. 
(c)). It can be inferred, from the agreement with the 
experimental magnetic moments, that the present nu- 
cleon g factors are reasonable. At A = 18, the SDPF-M 
and USD calculations still give similar magnetic moments 
in good agreement with the experiment. On the other 
hand, the quadrupole moment by the SDPF-M is larger 
by about 30% than the USD value. Recently, a very 
precise measurement of the quadrupole moments for Na 
isotopes has been carried out by Keim et al. 112, 113 • The 
measured quadrupole moment of ^^Na is 8.6(3) e fm^ in a 
good agreement with the SDPF-M prediction, 9.1 efm^. 
On the other hand, the deviation of the USD result for 
^^Na from experiment seems somewhat larger than that 
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FIG. 1: (a) Electric quadrupole moments, (b) magnetic dipole 
moments, and (c) npnh (n — 0, 2, 4) probabilities of the 
ground states of neutron-rich Na isotopes, as a function of 
the neutron number, A^. In (a) and (b), the circles are ex- 
perimental values taken from ,22. 23] , while the solid and the 
dashed lines denote, respectively, the MCSM calculation with 
the SDPF-M interaction and USD-model calculation. 



for typical sd-shell nuclei [7| . The situation is almost un- 
changed if the radial wave function is replaced with the 
Hartree-Fock one or if the isovector effective charge is 
tuned "2^. It was thus suspected in 23] that the experi- 
mental quadrupole moment of ^^Na might indicate some 
influence from the intruder configurations. The present 
MCSM calculation indeed shows, in Fig. ^(c), the large 
mixing of intruder configurations by ~ 42%, and their 
effects are visible in the quadrupole moment. 

Unlike the cases for N < 19, in the cases of = 19 
and 20, the moments cannot be reproduced by the USD 
model at all. Let us start with the most unstable iso- 
tope, ^^Na. The ^^Na nucleus {N = 20) has been known 
as a typical case of the intruder dominance in the ground 
state [3 • Its magnetic moment was reproduced by previ- 
ous shell-model calculations in a large shell-model space 
by Fukunishi et al. [35l | and by Caurier et al. 0, sup- 
porting this picture. As shown in Fig.^Cc), the present 
calculation, allowing full configurations within the va- 



lence shell, confirms the intruder dominance in ^-'^Na and 
indicates some mixing of even higher intruder configura- 
tions. Accordingly, we can reproduce not only the mag- 
netic moment but also the quadrupole moment |23| . Note 
that the prediction of the energy of the first excited state 
of •^^Na is in agreement with the measurement by 
intermediate-energy Coulomb excitation p^ . 

We shall now move on to ■^"Na [N = 19), which is 
the most crucial nucleus in this paper. The ground-state 
property of ^°Na had been rather obscure so far. The ob- 
served two-neutron separation energy shows no deviation 
from the USD-model systematics. The ground-state spin 
J — 2 [s^ can be explained by the USD model. This 
is in contrast to the anomalous J = 3/2 ground state in 
■^^Na, which is not obtained by the USD model. The ex- 
perimental magnetic moment of '^'^Na 2.083(10) ^jln |33l| . 
however, deviates from the USD-model value 2.69 /iAf. 
This deviation seems to be somewhat larger than typ- 
ical deviations in sc?-shell nuclei. As Fig. ^ (b) shows, 
this deviation is resolved by the MCSM with the SDPF- 
M interaction as a consequence of the intruder ground 
state (see Fig. ^ (c)). As the calculated magnetic mo- 
ments of ^°'^^Na can be changed only less than by 0.1 
/iAT by replacing the effective g factors with the free nu- 
cleon ones, the agreement with the experiment should 
not be attributed to the choice of the g factors. Re- 
cently, the quadrupole moment has been measured also 
by Keim et al. ,22!|. This value, even its sign, turns out 
to be quite different from the USD prediction. Figure 
(a) indicates that the MCSM with the SDPF-M inter- 
action indeed reproduces this quadrupole moment, too. 
Therefore, the properties of the electromagnetic moments 
indicate that, in Na isotopes, the ground state is domi- 
nated by the intruder configurations at iV = 19 ('^°Na), 
and intruder configurations are substantially mixed in 
the ground state already at iV = 18. 

It may be of interest to discuss the binding energies of 
Na isotopes to some detail, because the USD model ex- 
plains the observed trend of binding energies up to iV=19 
rather well. Figure [21 compares the two-neutron separa- 
tion energies (5*2^) of Na isotopes between the experi- 
ment 1^ 13 and the shell-model calculations. It can 
be seen that the USD gives an overall agreement with 
experiment as well as the MCSM calculation with the 
SDPF-M interaction, except for the failure by USD at 
N = 20. This problem at = 20 has been known for 
many years, as discussed in Sec. In fact, as previous 
models (see, e.g., 0)13)0)0) indicated, the shortage of 
the S2n of ^^Na by 1.5 MeV in the USD model is reme- 
died by having the intruder-dominant ground state. On 
the other hand, the 5*2,1 value at TV = 19 can be repro- 
duced well by both the USD and SDPF-M, whereas their 
wave functions are completely different as can be seen in 
Fig. ^(c). We shall now resolve this puzzle of ^°Na. 

Figure (a) compares the experimental S2n of '^'^Na 
with the calculated values by using the USD interaction 
and the SDPF-M interaction. With the SDPF-M inter- 
action, we carry out two calculations, i.e., a truncated 
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FIG. 2; Two-neutron separation energies of Na isotopes, as 
a function of the neutron number, A''. The circles and the 
crosses are the experimental values taken from the mass table 
by Audi et al. |33| and a new measurement by Lunney et al. 
[21|. respectively. The solid line denotes the MCSM calcula- 
tion with the SDPF-M interaction, while the dashed line the 
USD-model calculation. 




w 

in 
W 



10 



USD 

(b) shell ■ 



04, 



sd full 
^ SDPF-M / 



Exp. 



C16) 



16 



shell model within the sd shell and the full calculation. 
The results from the USD and the SDPF-M within the 
sd shell show rather different 52,1 values of "^^Na, despite 
the same model space. In order to understand this differ- 
ence, the ESPE is considered as shown in Fig. O (b). In 
the SDPF-M interaction, the ESPE of the 0^3/2 for small 
Z is higher than that of the USD interaction. This dif- 
ference is a consequence of the shell evolution mentioned 
earlier, and is the largest near Z — 8, creating a new 
iV = 16 magic number. The neutron 6^3/2 orbit is low- 
ered as Z becomes larger, due to the strong spin-isospin 
dependence of the A'' A'' interaction At Z = 11 (Na), 
this 0d3/2 is still rather high. Thus, if the calculation is 
restricted to the sd shell, the SDPF-M interaction pro- 
duces smaller S2n than that of the USD for the nuclei 
where the last neutron is in the 0^3/2. On the other 
hand, the intruder configurations dominate the ground 
state in the full calculation by the MCSM, increasing the 
binding energy and making S2n larger to the same ex- 
tent as the USD calculation in the sd shell. Thus, almost 
the same separation energies can be obtained from differ- 
ent mechanisms, and one has to combine other physical 
observables to draw definite conclusions. 



IV. ENERGY LEVELS 

The energy levels of ^^^■^''Na are calculated for the 
SDPF-M interaction by the MCSM, and are compared 
with both the experiment and the USD model in Fig. 0] 
Note that those of '^^Na have been reported in (13], and 
are not included here. Although there have been just 
few experimental levels published so far, they provide us 
with important information. We shall present, with em- 



FIG. 3: (a) Sin of ^°Na compared among the shell-model 
calculations (with the USD interaction and the SDPF-M one) 
and experiment. For the SDPF-M interaction, a truncated 
calculation within the sd shell and the full one by the MCSM 
are compared, too. The circle and the cross are experimental 
data taken from [3^ and [23], respectively, (b) Corresponding 
dominant neutron configurations of the ground state and the 
ESPE's obtained from each interaction. All the ESPE's are 
obtained by assuming the filling configuration. 



phasis on the intruder configurations, predictions from 
the SDPF-M interaction, which can be some help for fu- 
ture experiments. 



1. ^^Na 

The USD model and the MCSM with SDPF-M give 
similar energy levels, and in the latter a state substan- 
tially affected by the intruder configurations does not 
appear low. The calculated energy levels are in good 
agreement with the experimental ones observed recently 
[33 except for the absence of the 1.725 MeV state (see 
Fig. The agreement with the experiment confirms 
high predictive power of the USD interaction for near- 
stable nuclei. The state absent in the calculations has 
been tentatively assigned as the l/2~ 36]. As discussed 
in fsol], it is unlikely that at iV = 16 the negative-parity 
state dominated by a one-neutron excitation across the 
iV = 20 gap appears low, partly because of a somewhat 
large gap from the lsi/2 to the above orbits and partly 
because of the strong pairing correlation in even-A'^ neu- 
trons. If this state has a negative parity, it would involve 
a one-proton excitation from the Z = 8 closed shell. 
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FIG. 4: Comparison of the energy levels of ^^~^°Na relative to the experimental ground state among the experiment (Exp.) and 
the shell-model calculations by the SDPF-M and the USD interactions. The E2 strength from the ground state is illustrate d by 
the width of the arrow. The experimental B{E2) values of '^^'^'^Na and the energy levels of '^'^Na are taken from |25t and |36(I . 
respectively. For ^°Na, the levels calculated from SDPF-M interaction are grouped into four columns; the first (second) one is 
K=2 (1) rotational band dominated by intruder configurations, the third one represents spherical states which are basically of 
normal configurations, and negative-parity states are shown in the fourth column. 



2. ^^Na 

The experimental ground state of ^*Na is J = 1 [33l |. 
while in the calculations the 1^ and 2^ states are located 
quite closely in energy but the 2^ is slightly lower. The 
1^, 2^, 3^ and 4^ states are dominated by the configura- 
tions consisting of a neutron 2^(0^3/2)^ coupled weakly to 
the proton J = 3/2+ or 5/2+ state (see the energy levels 
of ^^Na), and are close to one another. For these states, 
both the shell-model calculations give similar excitation 
energies. 

A recent Coulomb-excitation experiment by Prity- 
chenko et al. shows a 7 ray at 1.24 MeV with 

B{E2) t= 54(26) e^fm^. The MCSM with SDPF-M 
gives b{e2;\+ 2+) = 69 fm^ and B{E2;l+ 
3^) = 47 fm^, either (or the sum) of which may cor- 
respond to the observed 7 ray. On the other hand, the 
B{E2; 1+ ^ 2+) and B{E2; 1+ 3f) are as small as 
19 and 27 e^fm**, respectively. Similar B{E2) values are 
obtained by the USD model, but the relevant 2^ and 
3^ energy levels by the USD are lower by 0.3 MeV 
than those of the SDPF-M (see Fig. By analyzing 
the occupation numbers of the wave functions, it turns 
out that the 2^ and 3^ states are mainly composed of 
one-neutron excitation from the lsi/2 to the 0^3/2- As 
the gap between these orbits is larger for the SDPF-M 



interaction, those states are pushed up. 

The negative-parity states are predicted to lie rather 
low reflecting the narrower TV = 20 shell gap, but there 
is no experimental information presently. They might 
be compared qualitatively to the state at 1.095 MeV in 
the = 17 isotone ^^Mg which can be a negative-parity 
state as discussed by Baumann et al. ^37] . 



The ground state of ^^Na is J = 3/2 experimentally 
[33I I. The calculations show very close 3/2+ and 5/2+ lev- 
els, and the MCSM gives the correct spin order, whereas 
the USD model does not (see Fig. This difference is 
because the 3/2+ state contains a larger mixing of the 
intruder configurations than the 5/2+. The shell-model 
calculations show that the 5/2^ state is strongly con- 
nected to the ground state with B{E2]2,/2"[ 5/2^^) = 
111 e2 fm* by the USD model (135 fm^ by the MCSM), 
while the B{E2) values from the ground state to the 
other normal-dominant low-lying states are very small. 
We thus point out that the Coulomb-excitation would 
hardly populate other excited states as far as the low- 
lying states are dominated by normal configurations. 

In the USD model, it is predicted that there are just 



6 



1/2^ and 9/2^ levels around 2 MeV. Apart from these 
normal-dominant states, the MCSM predicts, around the 
same energy, low-lying 3/2^, 5/2^ and 7/2^ states dom- 
inated by the intruder configurations. Due to the large 
mixing in the ground state, the 7/2^ may be excited by 
the Coulomb excitation with a moderately large value, 
B{E2;3/2f 7/2f) = 57 e^fm^, as predicted by the 
MCSM. 



4. 3"iVa 

Both the calculations succeed in reproducing the 
ground-state spin, but the energy levels are quite dif- 
ferent. In the USD model, the low- lying states are com- 
posed mainly of the configurations with a neutron hole 
1^(0^3/2)""^ coupled weakly to the proton J — 3/2 or 
5/2 state. The E2 strength between them should then 
be weak as depicted in Fig. ^ On the other hand, 
the MCSM with SDPF-M gives the intruder-dominant 
ground state which is strongly deformed. Indeed, we ob- 
tain a rotational band connected by strong E2 transi- 
tions: the E2 matrix elements calculated by the MCSM 
linked to the ground state are B{E2;2^ 3^) = 168 
e^fm^, B{E2;2+ 4+) = 90 e^fm^, and Q{2+) = 16 
e fm^ . They give rise to similar intrinsic quadrupole mo- 
ments, i.e., Qq = 58, 65 and 56 efm^, respectively, by 
assuming K — 2. The strong E2 transition has recently 
been measured by the Coulomb-excitation experiment by 
Pritychenko et al. ^5[: from the strength of the mea- 
sured 7 ray the B{E2 t) was deduced to be 130^^^ 

fm'' consistently with the MCSM calculation. The 
anomalous quadrupole moment (see Fig.^ and this large 
B{E2) value in ^°Na are excellently accounted for as a 
result of the large prolate deformation associated with 
the intruder-dominant configurations. 

From the viewpoint of the particle-rotor picture, the 
intrinsic state of the yrast band is regarded as a pro- 
ton in the 7r[211]3/2+ Nilsson orbit and a neutron in the 
i'[200]l/2+ coupled to a deformed ^*Ne rotor. As a re- 
sult, K = 1 and 2 are possible as the yrast band, and 
the MCSM shows that the latter is favored in energy. It 
is of interest to point out that this feature is consistent 
with the so-called Gallagher-Moszkowski rule [l^ that in 
strongly deformed nuclei the favored K is made so that 
the intrinsic spins of the last proton and neutron are par- 
allel. Thus, the agreement of the ground-state spin J = 2 
by the MCSM is not just an accidental fortune reflect- 
ing a particular interaction matrix element, but has been 
conducted once the intruder configurations dominate the 
ground state. 

The MCSM yields also the K = \ band starting at 
0.31 MeV. Its J = 2 and 3 members are calculated to 
lie around 1 MeV as shown in Fig. 01 while they are 
well mixed with the normal-dominant states. Also at 1 
~ 1.5 MeV excitation energy, normal-dominant spheri- 
cal states, corresponding to the lowest states in the USD 
model, appear as shown in Fig. 0] The negative-parity 



states are predicted to be rather low, dominated by the 
Iplh excitation across the iV = 20 shell gap. The com- 
petition between normal and intruder configurations in 
•^"^Na seems to be very intriguing, and is discussed in the 
next section in more detail. 



V. SHELL-GAP DEPENDENCE ON THE 
INTRUDER DOMINANCE 

From the above discussions on the moments and the 
levels, it becomes evident that the transition from the 
normal to intruder ground state occurs fully at iV = 19, 
after strong normal- intruder mixing already at N — 18. 
We shall show, in this section, that this normal-intruder 
transition for A^ < 20 is particularly sensitive to the shell 
gap. 

In general, an intruder state can be the ground state, if 
the energy gain due to dynamical correlations including 
deformation overcomes the energy loss in transcending 
nucleons across the shell gap. The shell gap is nothing 
but the difference between ESPE's of relevant orbits. The 
neutron ESPE changes rather gradually as a function of 
the neutron number, since the monopolc interaction for 
T = 1 is weak. Namely, the neutron shell gap is rather 
constant as a function of the neutron number. This im- 
plies that what is crucial in the transition from a normal 
to an intruder ground state within an isotope chain is 
primarily the neutron-number dependence of the corre- 
lation energy and its relative magnitude to the shell gap. 
Here, a good index of the correlation energy is the differ- 
ence between the eigenvalue of the total Hamiltonian and 
the expectation value of the monopolc interaction for the 
filling configuration. 

In Fig. 13 the sources of the correlation energy are 
sketched schematically. Since a normal state of a (neu- 
tron) semi-magic nucleus consists of configurations shown 
in Fig. |S1 (a), only the proton rearrangement is relevant 
to the correlation energy, which is generally small. On 
the other hand, the correlation energy is very large in the 
case of an intruder state composed of configurations like 
Fig. [51 (b), due to large numbers of particles and holes 
in active orbits. We note that the proton-neutron in- 
teraction produces much larger correlation energies than 
the interactions between like nucleons. This makes the 
correlation energy in Fig. [S] (b) much larger than that 
of Fig.|Sl(a), favoring the normal-intruder inversion even 
with a large shell gap. 

On the other hand, in the cases like Fig. [S] (c) and 
(d), a normal state of an open-shell nucleus has a neu- 
tron hole already. The neutron rearrangement is then 
possible, and strong proton-neutron two-body matrix el- 
ements contribute to the correlation energy. The intruder 
configurations of Fig. [3 (d) gain correlation energy simi- 
larly to the case of Fig|Sl (b). However, the difference of 
the correlation energy between Fig. [51 (a) and (b) is larger 
than that between Fig.[Sl(c) and (d), because of the sat- 
uration of the correlation energy with many particles and 



7 



normal intruder 

(a) (b) 



semi- 
magic 




(c) 





open- 
shell 



FIG. 5: Schematic sketch of the sources of the correlation 
energy of the intruder and the normal states of semi-magic 
((a) and (b)) and open-shell ((c) and (d)) nuclei. Typical 
configurations for these states are shown. The proton-neutron 
interaction is illustrated by thick wavy lines, while the proton- 
proton and neutron-neutron interactions are drawn by thin 
wavy lines. 



many holes as is the case in Fig. El (d) . A concrete exam- 
ple can be found with the USD interaction: a semi-magic 
■^'^Na gains the correlation energy only by 1.7 MeV within 
the sd shell, whereas it increases to 3.7 MeV for "^^Na and 
further to 7.2 MeV for ^^Na. The correlation energy of in- 
truder states increases more slowly due to the saturation 
as mentioned just above. This implies that the intruder 
dominance in A'^ < 20 nuclei becomes less favored as N 
goes down from 20. Hence, if the normal-intruder inver- 
sion still occurs, it should be due to a narrower shell gap. 
We shall present a more detailed account on this point 
now. 

The SDPF-M interaction indeed gives a narrow N = 
20 shell gap (~ 3 MeV) for Na isotopes. Note that it still 
reproduces the large gap 6 MeV) of '^''Ca, owing to 
its monopole property. We now demonstrate how such 
a narrow gap of Na isotopes plays a crucial role in the 
intruder dominance in '^'^Na, by means of a simulation 
based on the argument just above. Namely, we vary the 
shell gap from the value given by the SDPF-M interaction 
to larger values, to see what happens. This can be done 
by changing the monopole interaction between the Od^ /2 
and the 6^3/2 as 



SVr 



T=1,0 
Od, /2,0c 



Xx) = -0.3a;, -f 0.7x MeV, 



(1) 



where denotes the monopole interaction between i 
and j orbits with isospin coupled to T • The param- 
eter X is to control the ESPE: a; = represents the sit- 
uation with the SDPF-M interaction as a starting point. 




=i- 2' 







n — 'f — 


(b) 




// — 




• • 








J — // — 




3 f 

present shell gap (MeV) 

FIG. 6: (a) Quadrupole moment and (b) magnetic moment 
of ^"Na as a function of the A'' = 20 shell gap (controlled 
by a parameter x in Eq. (0). The experimental [23 . |23| 
and USD-model ones are denoted by the squares and crosses, 
respectively, (c) 2p2h probabilities in the (positive-parity) 
ground states of ^^~^^Na. (d) Energies of the 2p2h- and 
Ipl/i- dominant lowest states of '^"Na (denoted by "2p2h" 
and "iplh", respectively) measured from that of the OpOh- 
dominant lowest state. The range of the shell gap giving the 
npnh ground state (n = 0,1,2) is indicated by nex. Note 
that the corresponding shell gaps of ^^'^^Na are, respectively, 
smaller and larger by 0.24 MeV than the one of 30Na. 



A larger x means primarily a lower neutron 0^3/2 level, 
i.e., a wider A^ = 20 gap in ^^^Na. Other effects are minor 
in this nucleus. 

Figure presents the variation of the ground-state 
properties of Na isotopes as a function of the gap thus 
varied. Results are mainly about '^"Na unless otherwise 
specified. The SDPF-M interaction (at a; = 0) gives a 
narrow A = 20 shell gap, i.e., 3.3 MeV for '^°Na. As x is 
increased, the gap becomes wider, and the "'"Ca gap 6 
MeV) is given by a; '--^ 2.6. An intermediate value a; = 1 
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reproduces the gap of USD (~ 4.3 MeV), implying that 
the USD includes some fractional effects of the current 
shell evolution 0. 

The quadrupole moment and the magnetic moment 
in Fig. |S1 (a) and (b) are almost constant up to 
the gap ~ 4 MeV, as x is increased from 0. But, it 
jumps to values comparable to that of the USD model 
around the 4 MeV gap. In order to see how this rapid 
transition occurs, Fig. (c) shows the probability of 
the intruder configurations in the lowest positive-parity 
state for '^°Na. As expected from the change of the 
moments, the dominant component of the ground state 
moves rapidly from intruder to normal configurations at 
the shell gap ~ 4 MeV. 

It is interesting to compare this transition in ■^"Na with 
the ones of ^^'■^-'^Na shown in Fig.|H|(c). Compared to the 
pattern of '^"Na, notable differences are that (i) the shell 
gap causing the transition is larger in ^^Na (~ 5 MeV) 
and smaller in ^^Na (~ 3 MeV), and (ii) the transition 
takes place more slowly than '^°Na. The former is because 
of the difference of the correlation energies in the normal 
states of ^^~^^Na discussed already, and the latter is be- 
cause even-A'^ configurations are strongly connected with 
the pair-excited states via the pairing interaction. This 
is the reason why the intruder dominance in '^'^Na has a 
particular importance to clarify the shell structure of Na 
isotopes, and we now confirm that the narrower shell gap 
due to the shell evolution [l^ plays a crucial role. Note 
that at the gap of stable nuclei 6 MeV) the intruder 
dominance does not occur even in '^^Na. 

We finally discuss the competition of the dominant 
configurations in ^°Na including a negative-parity state. 
Figure|Hl(d) displays the energies of the 2p2h- and Iplh- 
dominant lowest states measured from the energy of the 
OpO/i-dominant state, as the gap is changed. In Fig. 01 
we can see what happens as x is increased from 0. The 
2f state is close to the 3j~ but stays lower consistently 
with experiment. If the gap is made larger, the ground 
state is switched to a negative-parity state around at 
3.5 MeV, and persists for a while. At a larger shell gap 
~ 4.5 MeV, a competition between a positive-parity state 
and a negative-parity one is encountered again, where the 
former is dominated by normal configurations. Finally, 
after this competition the normal-dominant ground state 
persists. The "island of inversion" picture (g seems to 
correspond to the gap near the second competition (i.e., 
around 4.5 to 5 MeV): with the weak-coupling approx- 
imation, the Ihuj and 2huj states were calculated to be 
located at 0.306 and 0.776 MeV above the normal one, 
respectively 0- 

The competition of normal- and abnormal-parity 
states in ^ONa {N = 19) can be compared to a famous 
example of the parity inversion in ^^Be (A^ = 7). Both 
are related to the narrow shell gap, but we point out a 



large difference between them: the latter is considered 
as the competition between the OpOh and Iplh states, 
corresponding to the second competition in the present 
paper. Thus, in the case of the A^ = 20 region a more 
drastic event occurs in spite of the normal-parity ground 
state of ■^"Na, reflecting a further narrowing of the shell 
gap. 



VI. SUMMARY 

In summary, we have investigated where the disappear- 
ance of the magic structure starts in the isotope chain of 
Na referring to its mechanism. It is suggested that ex- 
perimental electromagnetic moments, energy levels, and 
B{E2) values of ^^Na with A^ = 19 clearly indicate the 
dominance of the intruder configurations in its ground 
state, by combining with a shell-model calculation using 
the MCSM. The present result is in sharp contrast to a 
previous speculation based on the USD model from 
the viewpoint of the binding energy, where the disap- 
pearance was supposed to occur right at A^ = 20. The 
same conclusion as this speculation was drawn by the 
"island of inversion" model. The difference between the 
present calculation and the previous models is mainly in 
the behavior of the effective A^ = 20 shell gap for small 
Z: the gap is substantially narrow (about 3 MeV) for 
Na isotopes with the present SDPF-M interaction. Nev- 
ertheless, owing to the monopole part of the SDPF-M 
interaction, in particular, its spin-isospin dependent com- 
ponent, the well-known 6 MeV gap is restored for *°Ca, 
as an example of the shell evolution in stable and unsta- 
ble nuclei [l^- The validity of this argument has been 
confirmed in this study quite transparently with Na iso- 
topes below A^ = 20 where the intruder states are shown 
to be, most likely, unable to beat the normal states with- 
out a narrower shell gap. 
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